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Stochastic network description provide useful information about the hnk between glass transition 
temperature Tg and network connectivity. In multicomponent glasses, this permits to distinguish 
homogeneous compositions (random network) from inhomogeneous ones (local phase separation). 
The stochastic origin of the Gibbs-Di Marzio equation is predicted at low connectivity and the 
analytical expression of its parameter emerges naturally from the calculation. 
Pacs numbers: 61.20N-81.20P 



Most inorganic solids can be made amorphous by va- 
por deposition onto cold substrates. However, only a 
very few of inorganic melts can be supercooled by a wa- 
ter or air quench to yield bulk glasses which solidify 
at the glass transition temperature Tg. Oxides as vit- 
reous silica {Si02) and chalcogenides (e.g. GexSei^x) 
represent some of the best-known glass formers in na- 
ture. There have been numerous efforts to understand 
the nature of glass transition, and to relate Tg to some 
easily measurable quantities. Tanaka has proposed a re- 
lationship between Tg and the mean coordination num- 
ber [n|. Gibbs and Di Marzio have developed a second- 
order phase transition model and obtained an empiri- 
cal relationship (GDM equation) between the transition 
temperature and the density of cross-linking agents in- 
serted inside a system of molecular chains 0. However, 
there is still no universal relationship between Tg and the 
glass network connectivity, satisfied by allmost all kinds 
of glass formers. The first attempt of such a quantita- 
tive description has been given only very recently \^. 
We present in this Letter several important results con- 
cerning the glass transition temperature variation as a 
function of connectivity in network glasses. 

Mean-field estimates of Tg work only at low connectiv- 
ity, but they do not work at high connectivity because a 
stochastic description fails. 

In the chalcogen limit (low connectivity), an adapted 
version of the GDM equation for chalcogenides M can be 
obtained analytically from the stochastic description of 
the network. Consequently, the parameter f3 appearing 
in B can be computed in this limit from the coordina- 
tion number for any glass system, with simple and elegant 
sum rules. 

The combination of the GDM equation in Q and the 
topological calculation of the parameter j3 yields the cor- 
rect trends in the Tg variation up to r = 2.4, in agreement 
with Phillips' constraint theory ||^. Moreover, it sug- 
gests that usual curve-fitting with GDM equation hides 
the stiffness transition at f ~ 2.4, when realized over the 
whole concentration range because of the occurence of 
chemical ordering. 

The prediction of the model is parameter-free and it 
can be easily extended from binary to ternary, quater- 
nary and multicomponent network forming materials. 
We should stress at this point that the model is not in- 
tended to describe the physics of glass transition, but 
rather the connectivity dependence of the temperature 
of this transition. To this end, the chalcogenide systems 



are of particular interest for the application of the model 
because of the minor role played by kinetics, compared to 
the larger role of connectivity in determining the value of 
the glass transition temperature Q. For the reader's con- 
venience, let us first sketch the main ideas of the stochas- 
tic model. 

If the network of a certain glass system form a random 
network, then one should be able to treat statistically 
with equivalent fashion different states of structural de- 
scription. For example, in a binary B^Ai^x structure 
the mean probability p^ — l/2[2pAA + Pab] of finding 
an atom A among randomly distributed bonds A — A 
and A — B (with probability paa and pab) should be 
equal to the probability p^ of finding it among a random 
distribution of A and B atoms. Thus, p\ = p\- We have 
excluded the possibility of a B-B bond which occurs only 
in a modifier-rich glass structure. The bond probabili- 
ties pij (with (i,j)=(A,B)) should be proportional to the 
concentration x and (1 — x) (or p\ and p%)^ a statistical 
weight Wij related to the coordination numbers of the 
atoms A and B {rA and rs), and a Boltzmann factor 
involving both the glass transition temperature Tg and 
the bond energies Eij. The probabilities paa and pab 
can be written as follows: 



PAA^-^{l-xfe~^--'^-^^ 



PAB = 
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Z normalizes the bond probabilities. The statement: 
1 — X = p^ = p\ can be solved in terms of the con- 
centration X, because p\ is constructed with it: 



X = 
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Also, equation (j^) can be made parameter-free (i.e. with- 
out involving Ae — Eab — Eaa) by condidering the limit 
X — Q (pure chalcogen) where T'g(O) = Tq- One gets from 
d): Ae = Eab ~ Eaa = /cbTq In [r^/r a] Finally and 
most importantly, the relationship (0) can be cast in a 
more compact presentation by performing the derivative 
of (0) with respect to Tg in the limit a; = 0, and inserting 
the energy difference Ae. We obtain then the following 
parameter-free slope equation which is the central result 
S to be used in the present work: 



dTg_ 

dx 



n 



X = 0,Ta=To 
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As r is defined by: r = r^(l — a;) + tbX, one can obtain 
the derivative of Tg with respect to f by replacing in all 
equations @-@ a; by (f - r^)/(rB - r^): 



dr 



To 



r=rA,Tg=To 



[tb -?'A)ln 



(5) 



The value of the coordination number rg can be deter- 
mined in most of the situations by the 8— ./V rule, where N 
is the number of outer shell electrons of the considered 
atom [0 (and of course, in the forthcoming, ta = 2). 
From the second part of equation (||), it is possible to 
demonstrate the stochastic origin of the heuristic GDM 
equation at low connectivity (f = 2) and to show that 
the modified equation proposed by Varshneya and co- 
workers Q is the correct expression as long as this equa- 
tion remains linearly extrapolated. The GDM theory of 
glass transition (based on equilibrium principles) is in- 
tended to describe the variation of Tg of long polymer 
chains of equal length, with chain stiffness produced by 
cross- linking agents ||] . The Tg of the cross-linked glass 
is suggested to behave as: Tg = To(l — kX)~^ where 
X is the cross-link density and n a constant. From the 
construction of the theory, is easily conceivable that it 
might describe chalcogenide glasses as well, because the 
initial selenium network is also made of long polymeric 
chains, and the modifier atoms as Ge or As should play 
the cross-linking agents. To this end, the GDM equa- 
tion has been recently adapted with success by Varsh- 
neya and co-workers in order to describe the Tg trends 
in chalcogenides g|: Tg = To/(l - (3{f - 2))'^ (which 
we will denote by VGDM equation). The parameter /3 is 
obtained from a least-squares fitting of the experimental 
data. Now, if one performs the derivative with respect 
to f of the first-order Taylor expansion (linear) of the 
latter expression in the vicinity of f = 2, one obtains 
(3Tq. This can be compared with the right-hand side of 
equation (^ and yields an analytical expression for (3: 
f3^^ = {tb — 2) ln[^]. In a two-component chalcogenide 
glass (A,B), the constant /3 has a topological origin and 
can be easily computed from the coordination number of 
the modifier atom B. 

We shall prove that the factor appearing in the ex- 
pression of P has a universal character and can be ex- 
tended to M-component glass system, yielding the value 
of the parameter /? for any sytem, to be inserted in the 
VGDM equation. This result will still be obtained in the 
limit where the VGDM equation emerges naturally from 
a stochastic description of the network. 

Before, we shall consider a glass system made of three 
different kinds of atoms (say A, B and C, with respec- 
tive concentration \ — x — y, x and y), one of them being 
the chalcogenide atom of the chain-like initial structure 
(when a; = and y — 0). The coordination numbers of 
the involved atoms are va = 2, vb and re- The average 
coordination number is f = tbx -\- rcy + 2{1 — x — y) and: 
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We can still identify the derivative of the first-order 
Taylor expansion of the VGDM equation in the vicin- 
ity of f = 2 with the right-hand side of (^), where 
consequently a; = and y — Q and where the quan- 
tities dx/dTg and dy/dTg have the form presented in 
equation (|) (i.e. Infrs/rAj/To and \n[rc/rA]/To). By 
identification, this leads to the analytical expression of 
the parameter /? in a glass made of three components: 
/3-1 = (rs - 2)ln[rB/2] + {re - 2)ln[rc/2]. The ex- 
tension to multicomponent systems appears to be quite 
natural. f3 has the same sum rules as the resistance in a 
parallel circuit in electrokinetics, i.e. it is the sum of the 
1/(3 of each related two-component system AB, AC, etc. 
Finally, for a system made of M different kinds of atoms 
with coordination numbers r^, we just have to sum up 
the M — 1 contributions (r^ — 2) ln[Y] in order to obtain 
the theoretical value of /3^^: 



1 
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J2{n^2)ln 
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In this notation, ttlm is the coordination number of the 
chain atom, equal to 2. However, when applied to ternary 
systems, the stochastic description with three kinds of 
atoms A, B and C is able to predict the VGDM equa- 
tion up to f ~ 2.3 with the correct value of (3, satisfying 
(0), starting from a system which satisfies p^ = p^ and 

Comparison with experimental data: In order to mini- 
mize the influence of the preparation techniques, we have 
carefully selected data of i) glass systems prepared with 
the same heating/cooling rate and ii) glass systems with 
more than five different compositions ||9[. The initial 
value To has been averaged over a set of data found in 
the literature (Tq oi v — Se has been taken as 316 K, 
of u — S* as 245 K) . We have plotted in figure 1 several 
experimental data on binary, ternary and quaternary sys- 
tems with the parameter (3 computed from the involved 
coordination numbers r^ and the VGDM equation (solid 
lines). The dotted lines correspond to the slope equa- 
tion (||) and shows the stochastic origin of the VGDM 
equation at low connectivity (f ~ 2). For other systems, 
we have performed a least-squares fit of the parameter 
(3 from the VGDM equation, denoted as (3exp- Then, we 
have compared the results with the predicted parameter 
(3pr obtained from (0). In the system Si-As-Ge-Te the 
involved coordination numbers of the modifier atoms are 
r, = (4,3,4). Thus, (3pr = (3 In 2 -I- In 3)-^ = 0.31, in ex- 
cellent agreement with the fit j3exp = 0.30 |12]. Examples 
displayed in fig. la show also that the slope (@) yields 
accurate trends in the variation of Tg with respect to the 
average coordination number f at low connectivity. 



Prom fig. la, we can obviously see that the addition of 
a two-coordinated atom (as tellurium or sulphur, r^ = 2) 
in multicomponent glass systems does not affect the value 
of the parameter /3 (because of the rate ln[ri/2]). A qua- 
ternary system which involves a two-coordinated atom 
can therefore be considered as a ternary system. The si- 
multaneous use of the VGDM equation and the stochastic 
prediction of (^ can give the value of the glass transi- 
tion temperature of any composition, at least for f < 2.4 
which is the limit imposed by Maxwell rigidity M. We 
invite the reader to check on this basis that the Tg of 
GaioBioSgo is about 325 K. For greater values of f, 
one has to take into account intermediate range order ef- 
fects such as the existence of rings [g , Q or chemical 
ordering. However, the stochastic model suggests also 
that usual curve-fitting with the VGDM equation [Q , jlj] 
does not permit to distinguish between floppy and rigid 
regions (limited by f = 2.4), whereas the combination of 
VGDM and (0) show direct evidence of the stiffness tran- 
sition between these two regions. This is due to the fact 
that f3 is evaluated from the rendom network description 
and it fails around 2.4. To illustrate this observation, 
we have plotted in fig. lb for the P-Ge-Se compound 
the least-squares fit of the VGDM equation (dotted line) 
iPexp = 0.66, realized over the whole concentration range 
Jlql) and the VGDM equation, combined with the sum 
rule (^ for Ppr = 0.55 (solid Hue). The deviation occurs 
at f = 2.43, consistently with constraint counting argu- 
ments M. We have also plotted in the insert of fig. lb 
the value P^xp computed from the VGDM equation, as a 
function of the average coordination number, in In-Ge- 
Se systems fll] . This clearly shows the threshold around 
f = 2.4 between the values /3 = 0.55 (predicted from 
the model by the random network picture) and P = 0.72 
(predicted from the model by the occurence of chemi- 
cal ordering, as we shall see below). A global fit of the 
VGDM equation (dotted line) with a unique value of /3 
would not have given this information. 

Effect of chemical ordering: In most of the chalco- 
genide systems, chemical ordering occurs for f > 2.4, 
and therefore stochastic description fails, as shown in fig. 
1 However, the description is still useful, if the network 
can be thought as a set of compound clusters inside a 
random network. 

Again, let us consider a ternary glass system 
Ai-x-yBxCy with the corresponding coordination num- 
bers rA — 2, rs and re- If we assume that the general 
tendency of the glass is to form a demixed structure of A 
and B in stoichiometric proportions, then we can rewrite 

the system as: {Br^Ars)x/rACyAi_y^(^rA+rB)x/rA ■ For 
the glassy matrix, this defines an effective concentration 
of C atoms y'^^-f = -. — -, — ^ — r-, — and the average coor- 

" l—x(rA+rB)/rA ° 

dination number of the glassy matrix is given by: 



If we proceed as before, i.e. performing the derivative 
with respect to Tg and looking at the limit (cc, y ^ 0) in 
order to identify with the expansion of the VGDM equa- 
tion, we can see that the corresponding parameter f3 is 
defined as P~^ = {rc — rA) ln[^]. In other words, the pa- 
rameter /3 of a ternary system which displays chemical or- 
dering can be computed by considering only the remain- 
ing two-component glass. We have checked the validity 
of this rule on a set of germanium incorporated chalco- 
genides. We have considered all possible stoichiometric 
demixed structures at the tie-line composition (e.g. in 
GcxAsySei-x-y, the possible structures are GeSe2 and 
As2Se3). Most of the ternary III— IV — VI systems such 



as Sb — Ge 



/3e 



X glasses (X 
= 0.78 120|; X=Te, P.^p 



S, Pexp = 0.61 119|; X=Se, 
= 0.79 H) behave as a 



single binary IV — VI glass (as Ge — ^e) with parameter 
close to Ppr ~ 0.72 when r > 2.4. This is explained 
by the presence of Sb2X3 clusters {X — S, Se, Te) in- 
side the remaining random network of Se-Se and Ge- 
Se bonds. Also, data on Sn — Ge — Se systems show 
that Pexp — 0.68, when SnSe2 clusters are assumed, 
0.72 ^^. Coming back to the illustra- 
te compound (insert of fig. lb), the 



close to Ppr - 
tive In — Ge 



behavior of Tg versus f can be explained as follows. For 
f < 2.4, the structure can be described by a random net- 
work of Se — Se, Ge — Se, In — Se and Ge — In bonds. 
For this system, the VGDM equation and (R) describe 
the Tg(f) behavior, with f3pr = 0.55 is computed from 
(0) with ri = (3,4). For f > 2.4, the network struc- 
ture looses its random character because of the occurence 
of In2Se3 clusters |18]. Thus, the system behaves as a 
pseudo-binary Ge — Se system with the effective concen- 
tration (|) and Ppr = 0.72. 
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FIG. 1. a) Comparison of the VGDM equation with the 
stochastic calculation of f3 (M) with experimental measure- 
ments in binary, ternary and quaternary glasses (solid line). 
The slope equation (El) corresponds to the dotted lines. Sys- 
tems Ge-Te- {Se), n = (4,2) and Ge - {Se), r, = 4 
[13]: /3pr = 0.72. System Sh-Ge-Te~ (Se), n = (3, 4, 2) 
and Sb - Ge - (Se), n = (3,4) [4]: fJpr = 0.55. System 
Ge- Br - (S), n = (4,1) [14]: P^r = 0.48. b) Occurence 
of chemical ordering in ternary selenides (data from [4], [16], 
[17]). The stochastic description yields for all Ppr = 0.55, 
since r; = (3, 4) in the three systems. The solid line rep- 
resents the combination of this description with the VGDM 
equation. The dotted line correspond to the VGDM equation 
with P fitted from the whole concentration range, as realized 
in [4]. For f > 2.4, the systems behave as a binary glass with 
Ti = 4 and Ppr — 0.72. The insert shows P^xp versus the av- 
erage coordination number f for the In — Ge — Se compound. 
All data sets have been displaced by 100 K for a clearer pre- 
sentation. 



